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In this paper we provide a unified framework for quasi-species evolution and stochastic quanti- 
zation. We map the biological evolution described by the quasi-species equation to the stochastic 
dynamics of an ensemble of particles undergoing a creation-annihilation process. We show that 
this mapping identifies a natural decomposition of the probability that an individual has a cer- 
tain genotype into eigenfunctions of the evolutionary operator. This alternative approach to study 
the quasi-species equation allows for a generalization of the Fisher-theorem equivalent to the Price 
equation. According to this relation the average fitness of an asexual population changes with time 
proportional to the variance of the eigenvalues of the evolutionary operator. Moreover, from the 
present alternative formulation of stochastic quantization a novel scenario emerges to be compared 
with existing approaches. The evolution of an ensemble of particles undergoing diffusion and a 
creation-annihilation process is parametrized by a variable /3 that we call the inverse temperature of 
the stochastic dynamics. We find that the evolution equation at high temperatures is simply related 
to the Schrodinger equation but at low temperature it strongly deviates from it. In the presence of 
additional noise in scattering processes between the particles, the evolution reaches a steady state 
described by the Bose-Einstein statistics. 

PACS numbers: 89.75.-k, 87.23. Kg,05. 10. Gg 



I. INTRODUCTION 

The intriguing relation between evolutionary dynamics 
and statistical mechanics [J-U has attracted the inter- 
est of classical evolutionary theorists. Fisher 0, Q and 
Kimura Q have related their results to the second princi- 
ple of thermodynamics and to the theory of gases. Inter- 
estingly, the relation between evolutionary theory Q and 
quantum statistical mechanics is emerging from a series 
of independent works 04HI that show a class of phase 
transitions occurring in the evolution of haploid popula- 
tions and other evolving complex systems described by a 
Bose-Einstein condensation. In haploid populations, this 
transition is the quasi-species phase transition 0, [l6r - 
|20| in which a finite fraction of an asexual population 
ends up having the same genotype if the selective pres- 
sure is higher than a critical value and the mutation rate 
is smaller than a critical value. Moreover, in a recent 
paper it has been shown that a condensation transi- 
tion in the Bose-Einstein universality class occurs also in 
the evolution of diploid sexual populations in presence of 
cpistatic interactions. When this condensation occurs, a 
finite fraction of pairs of genetic loci in epistatic interac- 
tions is fixed. 

The deep relation between evolutionary theory and 
quantum mechanics formalisms extends also to the dy- 
namical description of biological evolution. The quantum 
spin-chain formalis m I22t and the Schddinger equation in 
imaginary time [23l425| have been shown to solve mod- 
els of asexual evolution. The role of path integrals in 
describing the temporal evolution of populations (26l - [28j 
has been highlighted recently. 

Biological evolution is essentially a stochastic process. 
The relation between quantum mechanics and stochas- 
tic dynamics has been deeply explored over the years 



[29(. In particular the existing theory of stochastic pro- 
cesses for diffusing particles is described by Langevin and 
Fokkcr-Planck equations. Stochastic quantization ap- 
proaches define stochastic diffusing processes whose prob- 
ability density converges to path integrals of some specific 
quantum system (3fj| . l3l| . Stochastic quantization has at- 
tracted large interest in the physics community and has 
been shown to allow for large numerical simulations of 
quantum systems [32j . In this paper we propose a uni- 
fied framework for quasi-species biological evolution and 
stochastic quantization. 

In the first part of the paper we describe the evolu- 
tion of an asexual population with overlapping genera- 
tions in the limit of large populations. We show that 
it is possible to describe the evolution of a population 
in terms of eigenfunctions of the evolutionary dynam- 
ics, and this description allows for a generalization of the 
Fisher theorem in presence of mutations. We show that 
the variation in time of the average fitness of the popu- 
lation is proportional to the variance of the eigenvalues 
of the evolutionary dynamics. Therefore the average fit- 
ness of the population asymptotically in time reaches the 
highest eigenvalue of the evolutionary dynamics. This 
relation is very general and applies as well in absence of 
mutations, for low mutation rate and for high mutation 
rate. Therefore it is not in general related to the local- 
ization of the eigenfunction with maximal fitness on the 
fitness landscape. Moreover, in the first part of the pa- 
per related to biological evolution we review the results 
of the Kingman model which are related to the emer- 
gence of quantum statistics in the equation of biological 
evolution. These results are related to the quasi-species 
description of evolutionary dynamics that neglects ge- 
netic drift, which has significant relevance in presence of 
small populations and small mutation rate fiN <C 1 or 
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in the cases in which the population is fluctuating and 
undergoes population bottlenecks. 

In the second part of the paper we consider an en- 
semble of elementary particles which do not only diffuse, 
like in existing approaches of stochastic quantization, but 
undergo also creation and annihilation processes mimick- 
ing the dynamics of biological evolution. The probability 
density P{x, t) that a generic particle is at position x at 
time t follows the quasi-species equation first proposed 
by Eigen [l6[ to describe biological evolution. This allows 
us to define a new dictionary for a mapping from the evo- 
lutionary theory to the stochastic dynamics of particles 
undergoing the creation-annihilation process. The en- 
semble of particles is subjected to a stochastic Gaussian 
noise and a potential energy. The stochastic noise can 
be mapped to the mutations of biological evolution and 
the potential energy can be mapped to a Fisher fitness 
function. Moreover the stochastic evolution strongly de- 
pends on the parameter ft, that we call the inverse tem- 
perature of the system, playing the equivalent role of se- 
lection pressure in the biological evolutionary mapping. 
The probability of a certain configuration of the ensem- 
ble of particles satisfies the quasi-species equation and 
can be decomposed into eigenfunctions of the evolution- 
ary operator with a discrete spectrum. For low values of 
P, /5 <C 1 (i.e., low selection pressure) the quasi-species 
equation is directly related to the Schrodinger equation, 
while for /3 > 1 (i.e., high selection pressure) the quasi- 
species equation strongly deviates from the Schrodinger 
equation. In order to solve in one case the quasi-species 
equation for all possible values of the inverse temperature 
we explicitly solve the case of an harmonic potential. The 
evolution turns out to be dissipative, with a relaxation of 
the configuration to the fundamental state. Moreover we 
investigate the role of genetic drift in this model, show- 
ing that thanks to the gaussian nature of the noise the 
systems always reaches a stationary state. In particular, 
asymptotically in time the spatial distribution of parti- 
cles fluctuates around the fundamental state predicted by 
the solution of the quasi-species equation. The smaller is 
the number of particles the stronger are the fluctuations 
aournd the mean. Moreover we include an additional 
noise describing a scattering process of the particles such 
that after each scattering process a particle takes a ran- 
dom position distributed according to a fixed probability 
distribution. When we introduce this scattering prob- 
ability the ensemble of particles is proven to follow a 
Bose-Einstein distribution. Interestingly, the relaxation 
occurring in absence of the scattering process obeys a 
dynamical equation which has its equivalent in the gen- 
eralization of the Fisher theorem found in the mean-field 
treatment of the biological evolution presented in the first 
part of this paper. In this mapping the ground state of 
the ensemble of particles corresponds to the state of max- 
imal reproductive rate of the asexual population. There- 
fore the relaxation of the ensemble of particles to the 
ground state corresponds to the relaxation to the state 
corresponding to the maximal eigenvalue of the evolution 



of an asexual population. Finally, in presence of the scat- 
tering process the Bose-Einstein distribution is emerging 
naturally with similarities to the Kingman model [7| for 
the evolution of asexual populations with infinite number 
of genetic loci. 



II. EVOLUTION OF ASEXUAL POPULATIONS 
A. The quasi-species equation 

The genome of an asexual organism is formed by a 
single copy of each chromosome. If we indicate by 
i = 1, . . . , N a genetic locus, a given genotype is deter- 
mined by the allelic states {a} = (oi, 02, . . . ,<7j, . . . o~n) 
at each genetic locus i. The allelic state Oi at each ge- 
netic locus i can take 4 values corresponding to adenine, 
thymine, cytosine, guanine, i.e. oi — 1,2,3,4. The pop- 
ulation evolves under the drive of selection that favors al- 
lelic configurations corresponding to higher reproduction 
rate, and mutations that increase the genetic variation 
in the population. We assume that the reproductive rate 
M^({cr}) also called Wright fitness of a genotype {cr} is 
given by 



W({a}) = 



-PU({a}) 



(1) 



where U({o~}) is the Fisher fitness and /3 is the selective 
pressure. If j3 = every genotype has the same repro- 
ductive rate. If /3 ^> 1 the difference in the reproductive 
rate of genotypes having different U({cr}) is strongly en- 
hanced. 

We assume that at each time there can be a birth or a 
death process. We assume that the birth process depends 
on the reproductive rate ^({c}) = exp[— (iU({a})\ and 
that the death process is a random drift. If we define 
the probability P({o~},t) that at time t an individual 
has genome {<r}, the dynamical equation of evolution of 
P({a},t) is given by the mean-field equation, valid for 
large populations or growing populations, i.e. 



dP({a},t) 
dt 



-f>u {{<*'}) 



p(W},t) 



M wiKi i — -p(M^) 

(2) 

where the operator M/ ct \i/ ct /t. is defined in (|2l)|) and 
QiMlWY) is g ivcn b y 

Q{{°}\W}) = \[[{i-^{^<) + ^\ (3) 



where /z is the mutation rate. The partition function Z t 
in Eq. (0 is given by 

Z > = E E Q(W}\W})e- m{r7 ' }) P(W},t). (4) 

We observe here that the equation ^ can be reduced 
to the well-known quasi-species equation 0, [l7l - l20| if we 
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make a change of variables t t' with dt/Z t = dt' . Let 
us now assume to know the solution of the eigenvalue 
problem 



M 



WI{T'} 



e-^ u ^\ n ({a'})} = X n 7r n ({a}) (5) 



with X n describing the discrete spectrum of this prob- 
lem and TT n ({x}) the normalized cigenfunctions. If we 
decompose the function P({a},t) on the basis 7rn({c}) 
of eigenf unctions , i.e. 



(0) 



and if we make the self-consistent assumption that we 
know the function Z t , the dynamical solution of Eq. 
for the coefficients c n (t) is given by 



t (t) = exp [X n G(t) - 1] c„(0). 



(7) 



In Eq. (|7|) the function G(t) is defined through the func- 
tion Z t according to the equation 



G(t) 



dt'— . 

Zf 



(8) 



Using the definition for Zt given by Eq. (|4]) together with 
Eqs. © and ([7]), we can close the self-consistent equa- 
tions and uniquely determine the evolutionary dynamics 
of the population. Therefore the partition function Z t is 
given by 



(9) 



where the average (•) is performed over the functions 
c n {t) given by (0. Finally, using ([7]) we can derive the 
equation obeyed by the partition function Z t = (A), 



1 d{\Y 



2 dt 



= (A 2 ) - (A)' 



(10) 



This equati on g eneralizes the Fisher theorem of natural 
selection [J |27| ] in presence of mutations and describes 
the fact that evolution is an off-equilibrium process. In 
fact, (A) = (W({a})) is the average reproductive rate of 
the population. Eq. (fl"U|) expresses the fact that this av- 
erage reproductive rate, asymptotically in time, evolves 
toward the fundamental state Ao, as long as the envi- 
ronment does not change, i.e. the Fisher fitness func- 
tion U({a}) and the selective pressure remain constant 
in time. The rate at which the average fitness change is 
equal to the variance of the eigenvalues A„ over the dis- 
tribution c n (t). In the Fisher theorem, valid in absence of 
mutations, the average fitness increases proportionally to 
the variance of the individual fitnesses. Therefore we see 
that the generalization of the Fisher theorem [l[, equiv- 
alent to the Price equation [f|, can be done substituting 
the variance of the individual fitnesses with the variance 
of the eigenvalues of the evolutionary dynamics. We note 



here that since c n (t) sum up to one, but they are not nec- 
essarily positive definite, in presence of mutations, the 
fitness of the population can either increase or decrease. 
In any case, asymptotically in time the population is de- 
scribed by the eigenfunction associated to the maximal 
eigenvalue Ao. We observe here that the Price equation 
corresponds to the decomposition of the Eq. (|10[) in a 
term independent of the mutation rate and a term de- 
pendent on the mutation rate. Moreover we observe that 
the Fisher theorem is not only valid in presence of con- 
stant fitness landscape, considered in this paper, but is 
also valid in presence of fluctuating environments. In- 
deed most of the recent literature on the Fisher theorem 
relates to fluctuating environments, where t/({er}) is not 
constant. For recent results in these interesting aspects 
of biological evolution we refer the reader to the recent 
papers glUl]. 



B. The Bose-Einstein condensation in the 
Kingman model 

The Kingman model 0, Q is one of the most interest- 
ing stylized models of asexual evolution where the quasi- 
species transition is observed. In the framework of this 
model the quasi-specics transition can be exactly mapped 
to the Bose-Einstein condensation in a Bosc gas. One in- 
teresting aspect of the Kingman model is that the evolu- 
tionary dynamics reaches an equilibrium due to the con- 
stant drive of random mutations. In the Kingman model 
each individual is assigned a single real parameter e > 
determining its reproductive rate W(e), i.e. 



W(e) 



-fSe 



(11) 



Moreover, in this model, after each duplication, a mu- 
tation occurs with probability a and a new offspring is 
generated with random fitness e drawn from a given dis- 
tribution p(e). 

Therefore instead of writing Eq. @ for the distribu- 
tion P({a},t) that an individual of the population has 
genotype {a} at time t, we can write the equation for the 
probability density P(e, t\to) that a random individual in 
the population is associated with a given Fisher fitness 
e at time t and had the last mutation at time <o- The 
evolution of P(e,t\to) is given by 



dP(e,t\t c 
dt 



(!-«)■ 



~^P(e,t\t ) 



P(e,t\t ). (12) 



The partition function Z t in Eq. (|12j) is given by 

Z t = J dt J dee-^P(e,t\t ). (13) 

Asymptotically in time we assume self-consistently that 
Z t — > Z. Therefore, in this limit the probability P(e, t\to) 
that an individual has fitness e at time t under the condi- 
tion that the last mutation happened at time to, is given 
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by the solution of (|T2]) . i.e. 

P(e,t\t ) = ap (e)e[ cxp( - /3e)(1 - a) *- 1 ] (t - t «) . (14) 

Finally, integrating over to, we can evaluate the probabil- 
ity P(e, t) that an individual at time t has a given Fisher 
fitness e independently of to, i.e. 

ft 

P(e,t)= d< ap(e)e [oxp( " /3e)(1 " Q) *" l]( *" to) - (15) 



Therefore the steady state solution for Pg(e) reached 
in the limit t ~ > oo is given by 



P* B {e) = ap{e) 



1 



1 



e 0(e-p. B ) — 1 



(16) 



wher 



Z/(l — a) and the probability that an in- 



dividual has fitness e is determined by the Bose-Einstein 
distribution. Finally, using the definition Eq. (|13[) we 
can find the self-consistent equation that the constant 
Z = (1 — o)e~"^ B needs to satisfy, i.e. the normalization 
condition 



dep(e)- 



,p{e-HFs) _ 1 ' 



(17) 



If pie) vanishes for e — > 0, the Bose-Einstein integral in 
Eq. (fT?)) can be limited from above. As a result, at high 
enough selection pressure and low enough mutation rate 
the system might undergo a condensation phase transi- 
tion in the Bose-Einstein universality class. Below this 
phase transition a finite fraction of the individuals in the 
population shares the same genotype corresponding to 
the maximal fitness. This is one of the principal exam- 
ples that show the so called quasi-species transition: for 
low mutation rate and high selection a finite fraction of 
the population is found to have the same genotype. Be- 
low this phase transition the system is not stationary 
anymore and the population average fitness is strongly 
dependent on the statistic of records related to the oc- 
currence of largest fitness. Interestingly a similar phase 
transition occurs also in evolving ecologies [l4|, where 
the invasive species might strongly reduce the biodiver- 
sity, and in evolving models of complex networks [H, [l2| , 
where there might be the emergence of super-hubs like 
Google in the World- Wide- Web. 



C. The Fermi-Dirac distribution in presence of 
negative selection 

The Fermi-Dirac distribution can be obtained in pres- 
ence of negative selection by a similar mechanism that 
generates Bose-Einstein distribution in the Kingman 
model. This model is mostly interesting because of the 
underlying symmetry between Fermi-Dirac and Bose- 
Einstein distributions. We assume that each individual 
is assigned a parameter e describing its adaptability to 



the environment. The birth rate is one and new individ- 
uals are generated at each unit of time t with parameter 
e drawn from a given distribution p{e). The death rate is 
given by a random drift with probability a and by a neg- 
ative selection proportional to exp[/3e] with probability 
1 — a. Therefore the dynamical evolution of the proba- 
bility density P(e, t\to) that in the population there is an 
individual with parameter e born at time to is given by 



dP{e,t\tp) 
dt 



P{e,t\to) 



The partition function Z t in Eq. (|18[) is given by 



Z t = J dt Q j dee^P(e,t\t 



(18) 



(19) 



Assuming in the asymptotic limit of large times Z t — > Z 
to solve Eq. (fT5]) we find that the probability density that 
an individual has fitness e at time t under the condition 
that it was born at time to is given by 

P(e,t\t ) = p^eMMCi-aJiHCt-to). (20) 

The probability density that an individual at time t has 
Fisher fitness e independently of to is given by 



P(e,t) 



(iiop(e)e[ cxp(_,3e)(1 " Q) *" Q ] ( *" to) . (21) 



This systems reaches the steady state in the limit t — > 
oo with the equilibrium distribution given by the Fermi- 
Dirac distribution Pp(e), and in particular we have 



(22) 



with e^^ F = Za/(1 — a). This distribution is the Fermi- 
Dirac distribution. The self consistent argument is closed 
by determining self-consistently the equations that Z has 
to satisfy, i.e. the normalization condition 



dep(e) 



1 



e P(e-fj, F ) _|_ I ' 



(23) 



The duality between Fermi-Dirac and Bose-Einstein dis- 
tribution has been also recognized in the framework of 
evolving network models [l(| [HI and of models for evolv- 
ing ecologies flil ]. 



III. EVOLUTION OF AN ENSEMBLE OF 
PARTICLES 

A. Evolution of a large ensemble of particles 
described by the quasi-species equation 

We now study an ensemble of particles located in a one- 
dimensional space that undergoes a creation-annihilation 
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process. We map the evolution of the particles to biolog- 
ical evolution. In this mapping, the fitness function cor- 
responds to the energy of a particle, and mutations cor- 
respond to a stochastic noise. Finally we assume that the 
probability to find a particle at a certain space-time point 
for large ensemble of particles obeys the quasi-species 
equation. We simply replace in Eq. ^ the biological 
population by a large ensemble of particles and the indi- 
vidual genome by a point in continuous one-dimensional 
space. 

In the following, we define P(x,t) as the probability 
density that a particle is at position x at time t. We 
assume that this distribution, in the limit of large number 
of particles, follows the stochastic equation inspired by 
Eq. @ for biological evolution 



dP(x,t) 



M. 



dt Z t 
where the partition function Z t is given by 



-P(x,t) (24) 



Z t = / dx' / dxQ{x,x')e- 0U(x,)) P{x') (25) 



and the operator M^*, applied to a function 

e -PU(x)f( x ), ac t s as 



W*">f(a/j\ = / dx'Q(x,x')e-P u ^f(x') 



(26) 

where Q{x, x') describes the stochastic noise playing the 
role of mutations for the evolution of the ensemble of 
particles. The inverse temperature /3, with /3 > 0, plays 
in this stochastic dynamics the same role as the selective 
pressure in the biological evolutionary dynamics. For 
simplicity we assume that the position x at time t + dt is 
related to the position x' at time t by x = x' + rj where r\ 
is a noise with Gaussian distribution. Therefore we take 



Q[x,x') = 



1 



2ir/3D 



drj 5(x' -x-7])e 2D ? 11 (27) 



where Df3 is the variance of the noise. 

Introducing the Fourier representation for the delta 
function we get 



Therefore for (3 <C 1 the evolution Eq. ([24]) can be writ- 
ten as 



dP(x,t) 
di 



D d 2 P(x,t) 
2 dx 1 



+ U(x)P(x,t) 



(29) 



This equation is particularly interesting since it is linear 
in the Hamiltonian H(x, k) and can be put in relation 
with the Schrodingcr equation in imaginary time. In the 
next steps of our calculation we consider this particular 
limit (/3 <C 1) to compare the results of our evolutionary 
dynamics of particles with the solution of the Schrodingcr 
equation in quantum mechanics. 

In order to derive a specific solution of the stochastic 
dynamics described by Eq. (|24]l we consider the potential 
energy U (x) of the harmonic oscillator, i.e. 



U(x) = —mui 2 x 2 . 



(30) 



Assuming that time and space dependence of the proba- 
bility density P(x, t) can be factorized, we solve the fol- 
lowing eigenvalue equation for the space-dependent part 
of this probability, i.e. 



M, 



e-^VnOr')! = A„7r n (x) 



The normalized eigenfunctions are given by 
n n (x) = B„(x)e"ro^ 2 



(31) 



(32) 



where B n (x) indicates the polynomials of order n that 
solve this eigenvalue problem. Two possible values for 
e = ei and e = 62 solve this eigenvalue problem, i.e. 



£ 1,2 



(3Dmu> 2 



±Wl 



D^muj 2 



- 1 



(33) 



As the square root is always larger than one, one has t\ > 
and £2 < 0. The eigenfunctions are only normalizable 
for positive e — e\. Then one obtains the eigenvalues 



1 



(1 - ( 0e 2 )"+ 1 /2 



(34) 



Mr, 



-PU{x' 



/(*') 



= Af / dx' / dk 



Since the eigenvalues A„ are always positive we define 
the energy spectrum of the system as the set of values 
m(x,k)]+iH*-*')f( x >)E n given by 



(28) 



A, 



-0E„ 



where Af 



j_ 

271 



(35) 



2l ^p D and the Hamiltonian H (x, k) of 

this system is identified as H(x, k) = -jfc 2 + U{x). The 
action of the operator M x ,i on a function f(x) for j3 <C 1 
is given by 



Therefore, using Eq. (|34|) and Eq. (|35|) we can calculate 
the spectrum E n , given by 



-pU{x'] 



fix'} 



\ ) 



PE n = 

For p < 1, using Eqs. (|33]l . we find 



(36) 



Af dx' dk {1 - PH(x,k)} e k[x - xl) f{x') . 



(1 



^mu 2 D+0(l3) 



£2 



-Vrm^D+0{f3) . (37) 
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Finally, we can calculate the spectrum E n of Eq. (1361) . 
at the first order in /3, which is given by 



E n = VmZ? oj I n + 



1 



(38) 



This spectrum coincides with the spectrum of the quan- 
tum mechanical harmonic oscillator if we put D propor- 
tional to h 2 /m. Let us continue to study the stochastic 
evolution of every possible value of j3. If we decompose 
the function P(x,t) in the basis Tr n (x), 

P(x,t) =^c n (t)Tr n (x) , (39) 

n 

we obtain the self-consistent dynamical solution 



Cn(t) =atp[\ n G{t)-t}c n (0) 



(40) 



where the function G{t) is defined by the function Z t 
according to the equation 



G(t) 



d£ 



1 



(41) 



The solution of our evolutionary equation is given by 
the definition for Z t in Eq. (|2~5j) . the eigenvalues Eq. 
(pTTj) . and the structure of the solution of the probability 
density P(x, t) given by Eq. we find 



Z t =J2 X nC n (t) = (A) 



(42) 



where c n (t) are given by Eq. (|40 
dM]), dH]) it is easy to prove that 



Finally using Eqs. 



l d(X) 2 
2 dt 



((A 2 ) - (A) 2 ) 



(43) 



The partition function Z t = (A) therefore describes off- 
equilibrium dynamics. From the dynamical solution Eq. 
([40]) for the coefficient c n (t) it is easy to prove that the 
probability distribution P(x,t) converges asymptotically 
in time to the fundamental state 7i"o(a;) associated to the 
largest eigenvalue Aq. This phenomenon has a biological 
equivalent in Eq. (|10[) according to which a population of 
asexual individuals has a fitness that reaches the highest 
eigenvalue of the evolution operator. 

We observe that for the ensemble of particles the max- 
imal eigenvalue A = Ao corresponds to the ground state 
energy E according to Eq. (f36|) . Therefore, with time 
the particle system relaxes to this ground state. 

Finally, asymptotically in time, the average energy 
(U(x)) of the particles is given by the average of U(x) 
defined in Eq. (|30[) on the fundamental eigenfunction 
given by Eq. (f3"2"j) . Therefore we find 



(U(x)) 



1 muj 2 D 



(44) 



In the limit (3 <C 1 we find that the average energy (U(x)) 
is given by 



(U(x)) = -Twj 



(45) 



In the limit /3 > 1 instead, we find that (U{x)) is given 
by a classical type of expression, i.e. 



(U(x)) = i/3(M 2 



(46) 



where we have put D = h 2 jra. 



B. Stochastic effects 

The analytical treatment we have presented in the 
previous paragraph focus on the quasi-species equation. 
It is therefore a mean-field treatment that neglects the 
stochastic fluctuations present in the dynamics that we 
are proposing to study. In biological evolution these ef- 
fects are known as genetic drift. Typically these effects 
are important for small populations iV and small muta- 
tion rates /x, i.e. Nfi <C 1. For example in the Kingman 
model we might have that the fixation time of the fittest 
mutant can be smaller than the typical time needed for 
new mutations to occur. When this occurs the evolution- 
ary dynamics strongly depends on the statistics of records 
at which the fittest mutants arise in the population. 

In this section we study the stochastic effects in the 
model of evolution of the ensemble of particles under 
the influence of an harmonic potential U(x). In order 
to simulate the evolution of particle we adopt the Moran 
process @. We assume to have an ensemble of N parti- 
cles. Each particle i = 1, 2, . . . , N is assigned a position 
Xi. Starting from random initial conditions, at each time 
step we choose a particle i to generate a new particle. 
The particle i is choosen according to the probability 



n 



-PU{xi 



(47) 



The new particle will have a position x new chosen with 
probability 



Q{x , Xi ~) 



1 



e w 



2j3L> ( Xi X ) 



(48) 



Finally we choose a random particle j to be annihilated 
and we put 

Xj ^x new . (49) 

In Fig. [T] we show the asymptotic distribution for the 
ensemble of particles as a function of the noise ampli- 
tude D and the number of particles N in the ensem- 
ble. The fluctuations of these distributions are stronger 
for small ensemble of particles (N small) and high val- 
ues of D. Nevertheless always the quasi-species solution 
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FIG. 1: (Colour on-line) The asymptotic distribution of an 
ensemble of N particles undergoing a Moran process driven 
by an harmonic potential U(x) = ^muPx 2 with m = 1 and 
u = 1 and selective pressure f3 = 10. The distribution are 
collected after equilibration of the system for a time window 
of 400iV time steps. The solid lines indicate the theoretical 
expectations based on the quasi-species equation given by Eq. 
(I32|l . The smaller the population the larger are the stochastic 
fluctuations around the mean behavior described by the quasi- 
species equation. 



well captures the average of the distribution. In fact 
the asymptotic distribution fluctuates around the eigen- 
function ir (x) (given by Eq. ([52]) ) corresponding to the 
fundamental states of energy Eq. The suprisingly good 
agreement of the quasi-species approach reside on the 
fact that this stochastic process is defined over a con- 
tinuous variable x and is described by a Gaussian noise. 
Therefore the diffusion can be arbitrary small but is al- 
ways present as long as D ^ 0. This is at odds with 
genetic mutations that might not occur over a finite time 
frame if the mutation rate is small enough giving rise to 
a separation of time scales in the system. 



to is given by 

dP(x,t\tp) ^ M x , x ,e~^P(x',t\t ) 

j t = (1-a) ^ P{x,t\t ). 

(50) 

where M Ij3: / is given by Eq. (|26[) and Q(x, x') is given by 
Eq. (|27|). The partition function Z t in Eq. (f50|) is given 
by 



dt / dx' / dxQ(x,x')e- 0U ( x,) 'P(x' ',t\t ) 



(51) 

with Q(x,x') given by Eq. (|27[). In order to solve Eq. 
(|5"U|) we have to solve the eigenvalue Eq. (|5T|) for the in- 
teracting case. The eigenvalues and the eigenfunctions 
remain unchanged and are given by Eqs. (|36l) . (|32[) . If 
we decompose the function g(x) in the basis of eigenfunc- 
tions of Ma,,!' we get g(x) = J2 n 9n^n(x) . Asymptoti- 
cally in time, the dynamics will reach a stationary state 
determined by a constant limit value Z t — > Z. We as- 
sume self-consistcntly to know the value of Z. With this 
assumption, at large times t we solve the dynamical Eq. 
([501) to find the probability density P(x, t\to) that a par- 
ticle is at position x at time t under the condition that 
the last scattering happened at to, i.e. 

P(x,t\t ) = Y j agne^ il - a) i- 1]{t - to) Trn(x) . (52) 



Finally, integrating P(x,t\to) over to we get the proba- 
bility P(x,t) that a particle is at time t at position x 
independently of to, i.e. 

P{x,t)= f dt o y;aflne [A " (1 - a) *- 1](t -* o) 7r„(x). (53) 
Jo 



Using Eq. (|53[) and neglecting the terms vanishing in the 
limit t — > oo we obtain that the probability p(n) of an 
eigenstate n is expressed in terms of the Bose-Einstein 
distribution 



C. Evolution in presence of a scattering process 



P{n) = g n a 



1 



oP(E n -ii B ) - 1 



(54) 



In the following we study the change in the evolution 
of the particles in the presence of a scattering process 
occurring with probability a where at any time there is 
a creation process. In order to mimic this process we 
assume that when a particle undergoes a scattering pro- 
cess it takes a random position distributed according to 
a probability distribution g(x) which does not depend 
on time. In order to study the stochastic evolution of 
this interacting system of particles we describe the evo- 
lution of the ensemble of particles that have not scattered 
since time to- The evolution for the probability density 
P(x,t\to) that a particle is at position x at time t under 
the condition that the last scattering happened at time 



which takes the form of a Bose-Einstein distribution with 
e -0»B = z/(l - a). Finally the constant Z = (1 - 
a)e~P^ B is fixed by the normalization condition 



i 



(55) 



This demonstrates that in this context it is possible to 
explain the emergence of the Bose-Einstein distribution 
by purely dynamical considerations. In particular this 
mechanism has its parallel in the Kingman model Q for 
the evolution of asexual populations with infinite number 
of genetic loci. 



8 



IV. CONCLUSIONS 

In this paper we have proposed a unified framework 
to study biological quasi-spccies evolution and stochas- 
tic quantization. We have shown that the dynamics of 
biological evolution of non-overlapping generations and 
large populations is naturally expressed in terms of eigen- 
values and eigenfunctions of the evolutionary dynamics. 
This approach allows for a generalization of the Fisher 
theorem in presence of mutations. In this framework the 
average fitness of the population changes in time propor- 
tionally to the variance of the eigenvalues instead of the 
variance of the fitness of the individuals. 

Moreover we have shown that the quasi-species equa- 
tion describing the evolution of a large ensemble of parti- 
cles provides an alternative path to stochastic quantiza- 
tion with respect to existing approaches. The probability 
distribution of the ensemble of particles can be decom- 
posed into the eigenfunctions of the evolution operator. 
In the low temperature limit the quasi-species equation 
is related to the Schrodingcr equation. Therefore, in this 
limit the spectrum coincides with the spectrum of the 
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